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Abstract 

We study the oscillations and conversions of relativistic neutrinos propagating in matter of 
variable density using the wave packet formalism. We show how the oscillation and coherence 
lengths are modified in comparison with the case of oscillations in vacuum. Secondly, we demon- 
strate how the equation of motion for two neutrino flavors can be formally solved for almost 
arbitrary density profile. We calculate finally how the use of wave packets alters the nonadia- 
batic level crossing probabilities. For the most common physical environments the corrections 
due to the width of the wave packet do not lead to observable effects. 
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1 Introduction 



Recent experimental discoveries [0 seem to suggest that neutrino oscillations P, ^, || really exist. 
Since the observation of neutrino oscillation may provide valuable information on the basic properties 
of neutrinos, e.g. masses and mixing angles, it is important to know the underlying physics also on 
the conceptual level. 

It was pointed out e.g. in Refs. ^ and ^ that the standard quantum mechanical treatment of 
neutrino oscillations using plane waves p, is not completely satisfactory for many reasons. The 
wave packet approach 0,[|]-[|15] provides a more physical picture which is particularly adapted for 
describing phenomena localized in space and time. This formalism also elegantly accounts for the 
loss of coherence by the separation of wave packets. However, some authors (e.g. [l^, 0) have 



been skeptical about the use of wave packets, and Refs. ||T^ |T8[ conclude that the concept of wave 
packet is unnecessary for all the relevant physical cases. A bunch of other methods has also been 
discussed [|^ |T6|, 0, (to name a few). 

In this paper we will consider neutrino oscillations and other phenomena in matter. We have 
decided to use wave packets because the calculations of the kind presented here have never been 
carried out before. 

We will first focus on the equation of motion for neutrinos propagating in matter of variable 
density. Generally the equation is modified in matter due to coherent forward scatterings ^ 
and an exact solution can be found only for few special cases. Here we address one special case 
where the density of matter changes slowly enough, so that the situation is said to be adiabatic. Now 
the eigensolutions for the equation of motion are found trivially, for arbitrary number of relativistic 



neutrinos. To describe neutrino oscillations in matter, we apply the method of Ref. |n] in connection 
with these solutions. We can express the results formally using effective oscillation and coherence 
lengths which are not local quantities anymore. 

Generally neutrinos may propagate nonadiabatically, and solving the complete equation of motion 
even for two neutrino flavors is usually far from trivial (see e.g. |2^, ^ for two specific density 
profiles). In this paper we show that the solution for "arbitrary" density profile can be constructed 
by using infinite integral series. Some supplementary calculations which may be of formal interest are 
enclosed in Appendix B. Nonadiabaticity is also related to the fact that so-called level crossings (or 
hoppings) between matter "eigenstates" take place (p4|, ^ and references therein). Our calculations 



show that the level crossing probabilities are modified when neutrinos are described by wave packets. 

2 Neutrino propagation in matter in the adiabatic limit 
2.1 Equation of motion 

Barring the details of the production and detection mechanisms of neutrinos, we we will focus on 
the propagation. Traditionally the propagation of neutrinos is modeled by solving the relativistic 
Schrodinger equation with the effective Hamiltonian 

H = ^p^ + + V{x, t)^p + Kn^p'^ + V{x, t) , (1) 
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where p is the momentum operator, m is the neutrino mass matrix and V{x,t) is a semiclassical 
potential due to the presence of medium. In this work we assume that V{x) does not depend on 
time, and that it is a matrix diagonal in the weak interaction basis or flavor basis. This is well 
justified when considering neutrinos from the Sun, from supernovae, or on Earth. On the other 
hand, this assumption excludes neutrinos in the early universe that must be treated otherwise. 

Instead of decomposing the wave function in momentum space, we prefer to consider the Fourier 
transform in energy space, 

ij{x,t) = f dEe-'^'ij{x,E) . (2) 



/27r 

This leads to a more physical picture when the Hamilton operator is a function of x. Particularly, 
now the energy is a constant, while the momentum p is a function of x. Throughout this work we 
will assume the propagation to be a one-dimensional phenomenon. 

The resulting time independent Schrodinger equation can be written in the relativistic limit as 



tdMx, E)=[-E+'^ + V{x)\ ^(x, E) + hl±^^\ ^(^, E) 



(3) 



where we have also assumed that E . The right-hand side can be locally diagonalized. We call 
the respective eigenvectors as local matter eigenstates, and the respective eigenvalues can be written 
as 



Pa{x,E)=E 



2E 



(4) 



where we introduced the effective mass fia{x, E) that simplifies the notation. In this paper the Greek 
indices refer to the flavor basis and the Latin indices to the matter basis. In the adiabatic limit, or 
for negligible mixing, Eq. (H) can be written as 



id^ipaix,E) 



-E + ^^j^] M^, E) = -pa{x, E)iJa{x, E) 



(5) 



Outside the adiabatic limit the complete equation of motion includes also nondiagonal terms (compare 
to Eq. (11) of Ref. @). 

From now on we follow the treatment presented in Ref. |T^. Eq. (|^) can be solved easily, and 
one has 
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dE exp 



; / dx'paix', E) — iEt 
Jo 



(6) 



where Eq. (|^) was used. The initial value of x is put to zero since we consider neutrinos produced 
at the origin. 



2.2 Wave packet solutions 

As initial condition, we assume that the neutrino wave function is a linear combination of matter 
eigenstates, 

mE)=Y.Ca{E)MO,E). (7) 
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In the relativistic limit the coefficients Ca match the respective row in the mixing matrix, Ca(x, E) = 
U*a{x,E), but for less relativistic neutrinos correction factors related to energy conservation must 



be taken into account ||3^. We assume further that each component of the wave function ijja{0,E) 
has a Gaussian form 

(E-E^r 



^,(0,E) = (27r4p)-i/^exp 



4a2 



EP 



(8) 



where Ea is the average energy of the corresponding matter eigenstate and aEP is the energy width 
related to the production process fulfilling the uncertainty relation aEpcrtp = 1/2. 
In the adiabatic limit we obtain for the partial wave packet the solution 



il)a{x,t) = Ci dEexp 



I I dx'pa{x' , E) — iEt 





(E-Eg 
^^Ip 



(9) 



where Ci is a numerical factor. We assume that the wave packet is sufficiently narrow in energy 
space, aEP -C Ea- This helps us to integrate the above integral and we also avoid considering the 
negative energies. This assumption is justified e.g. for solar neutrinos |T^, One can now expand 
the momentum as 



Pa{x, E) ^ Pa{x, Ea) + - Ea) = Pa{x, Ea) + 



(10) 



where Va{x,E) is the group velocity of each wave packet. Note that within this framework the 
group velocity depends on x via the potential V{x), unlike in some alternative works. Writing 
aEP ~ Va{0)app = Va{0)/ {2axp), one has after a simple integration 



ipa{x,t) = Caexp 



px 

/ dx'pa{x',Ea) - iEat 

Jo 



'■^ dx' 



lo Va{x') 



'IV. 



where Va{0) is the group velocity at the origin and app{axp) is the momentum (spatial) width related 
to the production process. Eq. ([TI|) implies that the wave function of a flavor neutrino is given 
by 



\u^{x,t)) = C,J2U:aiO)^^P 



[ dx'paix', Ea) - iEat - 

Jo 



^^Ip 



dx' 

lo Va{x') 



Wa) , (12) 



where U{x) represents the effective mixing matrix and the states \i>a) are orthonormal. 

To get out relevant physics we need to know the size of the wave packet. This is very non-trivial 
and lots of rather confusing estimates have been presented. Here we assume that the width of the 
neutrino wave is mostly related to the spatial details of the production process. In principle it also 
depends on temporal properties, like the stability of the state producing the neutrino, but this is 
relevant only for very short-lived particles. It has been estimated (0, and relevant references 
therein) that a^p is of the order 10~^ m for solar neutrinos, a^p ~ 10^^^ m at the neutrino sphere of 
a supernova, and a^p ~ 10~^ m for reactor neutrinos, for example. 
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2.3 Observation of neutrino oscillation 



The quantum state of the neutrino is measured by an appropriate reaction. Assuming the respective 
process being i/^ + X — (something visible), we can relate the quantum mechanical uncertainty of 
the detection process to the quantum mechanical state of the particle X before the collision. Here 
we assume that the relevant wave functions are Gaussians, centered at a distance L from the origin 
(source), with a spatial width axo- Hence the detection can be described by 

vl{L) ( r dx' ^ 



\Vl3{X 



L))=C',J2u;aiL)exp 



dx'pa{x',Ea) 



L Va[X' 



(13) 



where Va{L) is the group velocity at L. Note that this is independent of time fl^ 

We emphasize that the spatial uncertainty, cr^D, arises for purely quantum- mechanical reasons. 
One may be apt to believe that it is at most of the order of atomic distances, i.e. axo ~ 10~^ m — 
^g-io Yess. In reality there are also other instrumental uncertainties related to the resolution or 
dimension of the detector. These can be accounted for by taking an average over the relevant length 
scale. 

The amplitude of the process z/^ z/^ is now simply 

Aap{L,T) = 



dx{i'fs{x — L)\i'a{x, T)) 
C'^T^KamM [dxexp 

a •' 

^2(0) / r dx' 



Va{x') 



T 



dx'paix',Ea) 

J 

vl{L)_ ( r dx' ■ 



lEaT 



(14) 



The first two terms in the exponential are independent of x, but the integration of the two other 
terms is far from obvious. Notice that the analytical form of Va{x) is unknown in general, and that 
we have not even chosen any specific density profile so far. It turns out, however, that the integral 
can be evaluated with a saddle point method to a sufficiently good approximation. 

After a straightforward, but rather lengthy calculation, presented in Appendix A, we obtain the 
amphtude (with Eqs. (^) and ( [ATO| )) 



dxpa{x, Ea) - iEoT - - 
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+ 



_£xD_\ " 

ML)) 



(15) 



The respective probability of the process Va is then 



\Aa(3{L,T)\' 



(16) 



a,b 



where 



G(L,r) 



dx {pa{x,Ea) 



T 



4Xa 



dx 



Vaix) 



Pb{x,Eb)) - i{Ea 

2 

I / 

T 



1 

4x1 



Ei,)T 
L dx 



vi,{x) 



(17) 
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with 



Xa — 



O'xP 



+ 



<7xD 

ML), 



We still have to perform an integration over time in Pai3{L,T) since we are mainly interested in 
PapiL) 14|. A straightforward Gaussian integration of the exponential of Eq. ([T6|) yields 



~ exp < z / dx{pa{x,Ea) - Pb{x,Eh))\ 



exp 



X exp 



^ dx 



^ dx 



XaXb 



Xa + Xb 

\ y« Jo 



— {Ea — Eb 



1 I 1 dx 1 dx 

+ 



4 \XaJO Va{x) Xb -Jo Vbix) 



1 dx 1 dx 
+ — 



Xa Jo Va{x) Xb Jo Vb{x) 



(19) 



without the prefactor, which is a constant in the relativistic limit. A detailed analysis renders Eq. 
(|l)to 

J rfTe^(^'^) ~ exp |z dx{p,{x, Ea) - Pb{x, Eb)) - i{Ea - Eb)L 



n 2 



dx{Va{x) - Vb{x)) 



(Eg — Eb 

8a2 



(20) 



where = a^p + a^p,, ap = l/(2cr^) and the relativistic limit is once again considered. 
In order to simplify Eq. (pO]), we approximate (compare to Ref. [Q) 



Ea^Eo + ^ 



2Ea 



(21) 



where /(/^(O) is the effective mass at the origin, Eq is the central energy in the limit of zero neutrino 
masses and ^ is a parameter of order unity, related to the energy-momentum conservation of the 
production process. The corresponding momentum is given by 



Pa{x, Ea) ^Ea- ^ ^0 + ' ^ 



2Ea 



2E, 







2E, 



(22) 







Notice that Ea is a constant, but Pa{x, Ea) depends on x, thus this approximation does not contradict 
the chosen convention. One may also point out that the initial dependence = fj.a{x, Ea) has reduced 
to fia = fJ^a{x,Eo). Using Eqs. and (0), Eq. (g^) gives 



J dTe^^^'^^ e^V I' dxAt^lix) - ^ 



at. \Jo 



dxAvab{x) 



{Eg - Eb)' 
Sal 



(23) 
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where A/i^^(x) = fii{x) — fJ^Kx), Avab{x) = Va{x) — Vb{x) and the last term has not been altered on 
purpose. 

Combining the above results we obtain a compact formula for the probability to observe a neutrino 
in a state at a distance L (with the correct normalization J2f3 PapiL) = 1), 



a,b 



X exp 



-2m- 



{Ea - E,f 



(24) 



'LZ{L) \L2\L)) Sal 
where the effective oscillation and coherence lengths are defined by 

Jo f^a;A/i^f,(x) dxAvab{x) 

The problem has reduced to computing the integrals over the effective mass and group velocity 
differences along the neutrino path. 

Let us comment on the following issues in connection with Eqs. (|^) and (pSf): 

1. The fact that the effective mixing matrices depend on the production and detection location is 
one part of the well-known MSW effect pi], 



2. Compared to the usual oscillation and coherence lengths, L°^^ = AttEq/ Aq, L^°'^ = 2\/2ax/\Avo\ 
(where Aq{Avq) is the mass squared (velocity) difference in vacuum), we see that the corre- 
sponding effective lengths take into account the changes of A/i^^(a:) and Avab{x) over the whole 
path of propagation. In this sense the oscillation and coherence lengths are not local quantities 
anymore. 

3. The physical coherence length, corresponding to a length scale where the oscillation ceases, can 
be obtained by solving the equation L = L'^l^{L). The physical oscillation length for variable 
density lacks a clear definition. 



4. The last term of the exponential in Eq. ( plf ) is related to the energy conservation within the 
uncertainty o"p. Its physical meaning is easy to understand: if e.g. \Ei — i?2| ^ o-p, only one of 
the states ui or 1/2 is "allowed" , i.e. there is no oscillation. 

5. The treatment presented here is not limited to any specific density profile if only Eq. (^ is 
valid (i.e. the adiabaticity is in effect, the matter density is not too high etc.). 



The discussion of Ref. [|1J] is in many respects applicable also here. 

We finally point out that the calculation of Jq dxAvab{x) is trivial in the relativistic limit: the 
group velocity is given by definition by 

dpa 2pl{x) + 2p„(x) dp^ ' ^^^^ 
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and taking into account the approximation of Eq. 



+ (27) 



Hence 

f-^ . . . . 1 



/ dxAvabix) ^ --^{-1 + EoOe,) f dxAfil^ix) , (28) 
JO 2i^o "'0 

i.e. we see that dxAvab{x) can be obtained easily if Jq dxAjj,l^{x) is known. 

3 Examples 

As an example we calculate the effective oscillation and coherence lengths for two specific density 
profiles. For simplicity only two neutrino flavors, Uf, and u^, are considered. 
The effective mass squared difference in matter is given by 

Amix) = \l (Ao COS 20 - 2V2EQGFNe{x)y + sin^ 20 , (29) 

where Aq is the mass squared difference in vacuum, 9 is the vacuum mixing angle, Gi;' is the Fermi 
constant and Ne{x) is the number density of electrons. (For other neutrino flavors or exotic matter 
contents Ne{x) should be replaced by an appropriate combination of the particle densities of the 
matter.) The resonance, where the effective mixing is maximal, is reached at the density 

An cos 26 , , 



The passage through the resonance is adiabatic when 

Aosin2 20 N,[xn) 



En COS 29 



» 1 , (31) 



where Q is the adiabaticity parameter. It is also assumed that \GFNf,{x)\ <C Eq. 

One should notice that the matter may affect considerably the oscillation of neutrinos propagating 
in constant density. This is the case e.g. if iVg is very close to its resonance value and 6 is small. 
Then (see Eq. (pQ])) A^ <^ Aq and the oscillation length may be highly longer than in the vacuum 
case. Since at the surface of Earth p ~ 3 g/cm^ and ~ 1/2, and hence GpNg ~ 10~^^ eV, the 
effect might be observable in long baseline experiments e.g. for Eq ~ 10 GeV, Aq ~ 10"'^ eV^ (and 
small 6). Also the group velocity difference, and consequently the coherence length, can be modified 
for suitable parameter values (cf. Eqs. (15) and (16) in Ref. |]TD[)- At the resonance the coherence 
length is usually increased. 
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3.1 Linear density profile 

The linear profile is by far the most important example, because many actual profiles can be locally 
approximated by it. Let us parameterize the density as 

iVe(x) = \{k - x) , (32) 

where A and k, are parameters. The adiabaticity condition, Eq. (pi]), leads to 

(Ao sin 26)^ 



We first calculate 



> 1 . (33) 



with (see e.g. |3 



^rfxA„(x) = - I{k - L)] , (34) 

4c 



4ac - 6^ 



I{x) = {2cx + b)Va + bx + cx"^ H \n{2J c{a + hx + cx'^) + 2cx + h) , (35) 

2a/ c 

and a = Ag, h = -4V2EoGi.AAo cos 261, c = ^{EQGF^f- A direct substitution to Eq. (ID then 
yields the effective oscillation length, in principle. Using Eqs. ([25|), (^) and (P^), the effective 
coherence length could also be calculated, but here that tedious calculation is omitted. 

Since Eq. ( |5^ ) is not very illustrative, we will next consider the low density limit (i.e. small A). 
One has 

/ dxA^ix) = LAo + V2EoGfX cos 26(1^ - 2kL) + Oi\^) (36) 

JO 



yielding 



■ V2EoGF\cos2e 
1 (L - 2k) 



(37) 



Ao 

where the second term in the brackets can be regarded as the first order correction due to matter. 
The coherence length is to this order (with Eqs. ( pSj) and (|2^) ) 



(38) 



i.e. the same as in vacuum. The lowest order correction of the coherence length is in fact always of the 
second order: the linear term ~ EoNe{x) (see Eq. (^)) is wiped out by the "operator" — 1 + EoOeq 
in Eq. (H). 

3.2 Exponential density profile 

Let the density of electrons be given by 

Ne{x) = Xe-""^ , (39) 
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where A and n are parameters. Eq. (^Tj) yields in this case 

Ao sin^ 26 



En cos 20k 



> 1 



Now we can write 32 



("^ dxA^{x) = -[I{l)-I{e--'^)] 

Jo K 



(40) 



(41) 



where 



Kx) 



Va + bx + cx"^ 



a In 



2a + bx + 2Ja{a + bx + cx"^) 



X 



+ —= ln(2y^c(a + 6x + ca;^) + 2cx + b) 
2v c 



(42) 



and a = Aq, b = —4:^/2EqGf^^oCos29, c = 8(-E'oG'fA)^. Once again, this result with Eqs. ( pSf ) and 
(|28|) would allow us to obtain the effective oscillation and coherence lengths in principle. 
In the low density limit we can expand 



dxAr^ix) 







K ' ' AqK 

Now we get the first order correction to the effective oscillation length 



(43) 







Ao 



1 - 



2V2EoGpX cos 26' 
LAqk 



-kL 



-1) 



(44) 



Eqs. (|25|) and (pSD yield the correction to the effective coherence length 



j^coh ^ 



|Ao| 



^ 2{EoGAsm29)\ .^.l _ 



(45) 



which is of the second order, in accordance with the remark made in the previous example. 

We conclude with a simple numerical application on solar neutrinos. The electron number density 



in the Sun is approximately [33 



X 



Ne{x) = 245A^Aexp y-^^-^^^ j cm" 



(46) 



where Na is Avogadro's number. Hence (Eq. (|39|) ) A ~ 10^^ eV^ and kL ~ 10^ if neutrinos are 
produced in the center of the Sun and detected on Earth. We use the values Aq ~ 10~^ eV^ and 
£"0^1 MeV, which fulfill the adiabaticity condition, Eq. (^0]), even for rather small values of 6. Eq. 
yields the correction of the effective oscillation length 



EoGfAcos20 



10~^cos2^ , 



(47) 
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where Gi? ~ 10 eV ^. The correction of the effective coherence length is similarly 

(^2^41^1^- 10- W2«, (48) 

where Eq. (|45| ) was used. It is thus seen that the matter effect modifies the oscillation and coherence 
lengths of solar neutrinos only slightly at least within the framework of this example. This fact 
is, of course, due to the insignificance of the density outside the Sun, i.e. the neutrinos propagate 
mainly in vacuum between the Sun and Earth. Similarly, it is to be believed that the matter effect 
is unimportant for the oscillation of all extraterrestrial neutrinos (excluding, of course, MSW effect 
and parametric resonance). 

We finally call your attention to the fact that even though both the linear and exponential density 
profiles violate the condition \GpNf,{x)\ -C Eq for some values of x, the sharpness of the wave packets 
forces X to be situated in harmless region (see Eqs. (|T2p and (|TBp). It is due to this same reason 
that the profile of Eq. (^61), written in spherical coordinates, is directly applicable. A more obvious 
reason is seen in the integrals in Eq. (pS]). 



4 Nonadiabatic neutrino propagation in matter 

4.1 Solution of the equation of motion for two neutrino flavors 

In the previous sections we have assumed that the neutrino propagation is adiabatic. In most physical 
environments, however, the matter density may change so rapidly that the nonadiabaticity must be 
taken into account. When neutrinos propagate nonadiabatically the local matter "eigenstates" are 
not anymore eigenstates but they become mixed. This means that transitions (i.e. level crossings) 
between the matter states may occur at certain probability. Mathematically this fact manifests itself 
in the complete equation of motion (in the matter basis), which includes also nondiagonal terms, 
whereas in the adiabatic limit the corresponding equation is diagonal (cf. Eq. (|)). 



In this section we look for a formal solution to the two fiavor equation of motion pj, ^ ^ 



'"^ [ M^, E) )-[ ^e'^{x) ^ ) [ M^, E) ) ' ^^"^ 

where A(a;) = iJi\{x) — fil{x) (this is not necessarily the same as in Eq. ( pOj ) if the intermediate 
matter contains e.g. muons), O'^ix) = dx6m{x) and Om{x) is the mixing angle in matter. In the limit 
4:E\9'^{x)\ -C A(x), equivalent to the previously considered adiabaticity condition, Eq. (^) (for 
^ oscillation), the nondiagonal terms can be neglected, and Eqs. (^) and (|^) (with a = 1,2) 
correspond to each other perfectly (discarding terms proportional to the identity matrix). Level 
crossings, on the other hand, are obviously caused by the nondiagonal terms ±i9'^{x). These issues, 
as well as exact calculations of the level crossing probabilities for specific density profiles, have been 
extensively discussed by numerous authors; see e.g. |pi]|,P^-[p5|],pi|-p5[ . 



One should remember that Eq. (^) describes the physical situation accurately enough if the 
following conditions hold (see e.g. Ref. 
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1. Neutrinos are relativistic. 

2. The density is low enough, i.e. \GfN{x)\ <C E, where N{x) is the number density of particle 
species relevant for a case under consideration (mentioned already in the previous section for 
N{x) = N,{x)). 

3. The density must not change appreciably over a length scale equal to the neutrino's de Broglie 
wavelength. 

The expression "arbitrary density profile" , used in the following, is to be understood in the context 
of the abovementioned limitations. 

We now demonstrate how Eq. (^) can be solved. Defining (we omit here the energy dependence 
of the wave functions) 

one has 

( ) = ( " , «W ) ( tf\ ] . (51) 

where B{x) = — iO'^ix) (and B*{x) = + iO'^ix))- The solution of Eq. (|5T|) is obtainable 

after some effort: 

M^) = CiH_(5,5*) + C2S+(i?,5*) , 

02(x) = C,E4B*,B)-C2E4B\B) , (52) 
where Ci and C2 are constants, and 

S±(y, Z) = 1±iJy-JyJztiJyJzJy + JyJzJyJz±--- (53) 

with e.g. JYJZJY^J:^ dx,Y{x,) J^^ dx^Zix^) dx,Y{x,). Eq. (0) yields finally 

i^iix.E) = ^Mxo,E)iE4B,B*) + E4B*,B)) 

-^Mxo,E)iE+iB,B*)-^+iB*,B)) , 

Mx,E) = ^Mxo,E){-E_{B,B*) + E_{B*,B)) 

+^Mxo,E){E+{B,B*) + E+{B*,B)) . (54) 

A straightforward substitution shows that this is really the solution of Eq. (^). As far as we know, 
no such complete solution has previously been presented in the literature. It is to be emphasized that 
Eq. ( p4D applies to arbitrary density profile, but on the other hand it has the undeniable deficiency of 
being somewhat formal and perhaps not very useful for practical calculations. Further discussion on 
the solution of Eq. (^) and on the corresponding equation in the flavor basis is found in Appendix 
B. 
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Let us show briefly that our solution gives a meaningful result in two specific cases. In the 

\Jx) 
AE 



adiabatic limit B{x) ~ ~ B*{x) and hence 



(55) 



and similarly 

iIj2{x,E) = ilj2{xo,E)exp ( -i [ dx^j^r] ; (56) 



AE 



i.e. the correct result is obtained. In the extremely nonadiabatic limil^, on the other hand, B{x) ~ 
—i9'^{x) ~ —B*{x) (in the resonance region), yielding 



^i(a;, E) = ^i(xo, E) (^1 + J B J B + ■ ■ ~ ^2{xo, E) (^i J B + i J B J B J B + ■ 
Mx,E) = ^^{x^,E){i j B + i j B j B j B + --^+^2{x^.E){l + j B j B + --^ . 
One may assume that 9m{x) changes abruptly from tx/2 to 9 (the vacuum mixing angle) in the 



B ^ -i(9 --] (5^ 



if Xq{x) is before (after) the resonance. Hence 



2 



TT 

sin ( - - 



i)i{x,E) = ?/'i(xo, -E) cos (^6* - - V^2(a;o,-£^) sh 
Mx,E) = Mxo,E)sm(^9-^^ +tlj2{xo,E)cos(9-^^ . (59) 
Putting e.g. ipi^xo, E) = 0, 4'2{xo, E) = 1, we see that the level crossing probabihty is 



\iJi{x,E)\^ = sm^(9--) =cos^9 , (60) 



as it should in this limit 



4.2 Level crossing probabilities and wave packets 

We now present how the existing results on the level crossing probabilities can be combined with the 
wave packet description in a consistent manner. As in the previous section, we restrict to two flavors 
which is sufficient for understanding the relevant phenomena. 

Consider a neutrino that propagates initially as tp2 and has the usual Gaussian form (Eq. (^), 

i.e. 

{E-E2f' 



ipi{xo,E) = and ip2{,XQ,E) = N exp 



Aalp 



(61) 



^Here we consider the usual textbook example where the matter states are related to electron and muon flavor 
states, and the matter does not contain muons. Remember that Eq. (^ ) is not necessarily limited to this standard 
case. 
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where = {^Inalpy^/^. From Eq. it follows that 



N 



exp 



{E - E2) 



4a2 



EP 



N 



i^2{x,E) = yexp 



{E — E2 



f{xo,x,E) 
h{xo,x,E) , 



(62) 



where f{xo,x,E) and h{xQ,x,E) represent the corresponding expansions in Eq. (|54D. Hence 



ipi{x,t) 



N 



2V27r 

where Eq. was used. Defining 



dE exp 



(i? — E2 



f{xo,x,E) 



(63) 



exp 



{E-E2)' 



4a2 



BP 



f{xo,x,E) = g{xo,x,E) 



one has 



?/'i(x,t) 



g{xo,x,t) 



(64) 



(65) 



where the circumfiex stands for the Fourier transform. The level crossing probability is \ipi{x,t)\'^ 
but we integrate over time (cf. Sec. 2) and have 



dt\g{xo,x,t)[-' 



dE\g{xo, X, E) 



(66) 



where Parseval's identity was used. We can express the level crossing probability for wave packets 
using the respective probability for a plane wave, Pic{E) = ^\f{xo,x,E)\'^. Hence 



Pic{E2,aEp) = |^i(x)|2 = / dEexp 



{E - E^f 
2a|p 



Pic{E) 



(67) 



where Pic{E2,(Jep) is the generalized level crossing probability that takes into account the energy 
width of the wave packet. 

The wave packet effects can be seen more clearly by expanding Pic{E) in series. Assuming ctep 
to be small (Sec. 2), it is sufficient to take the lowest terms, and Eq. ( |67|) gives (with E2 — > E) 



Pic{E, aEp) = Pic{E) + 



alp^POEl 
2 dE^ 



+ 0{^EP^^ 



(68) 



This equation clearly indicates that the use of the wave packets modifies the usual level crossing 
probabilities, Pic{E). It also shows that the wave packet correction is easily calculable if Pic{E) is 
known. 

Let us consider two simple examples: 
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1. Linear density profile. The well-known Landau-Zener probability is 



Plz{E) = exp ( --Q 



TC 



(69) 



where Q, given in Eq. (|3T|), should not be too small. Using Eq. ( |68D and remembering that 
Q ~ -p- (see Eq. (^)) one has 



(<yEpV 
\2E J 



3ttQ 



exp 



(70) 



For the most relevant cases (i.e. Plz{E) not too small) {j^QY /2 — SnQ is of the order of unity. 
2. Exponential density profile {Nf,{x) ~ e~'^^). Now j^S], ^ 

g-^5(i-cos2e) _ g-27r5 sinh(5 - A) _^ 



Pic{E) 



-2-K& 



sinh B 



where 



(5 



Aq 
2Ek 



A 



Tib 



1 - cos 20) , B = nS . 



A tedious calculation gives 



PUE,aEp) = Pic{E) 



1 + 



(71) 
(72) 

(73) 



where 

r = A{l + coth{B-A)){A-l + B{cothB-l))+B{coth{B-A)-cothB){l-BcothB) . (74) 

Again, for the interesting range of parameters F ~ 0(1) at most when e.g. solar neutrinos are 
considered (with E = 1 MeV, Aq = 10""^ eV^, and the density profile as given in Eq. (^)). 

The wave packet corrections turn out to be negligible for small aEp/E. This happens to be true in 
most physical environments, e.g. for solar neutrinos Gep/E ~ 10~^ — 10~^ [ [T^ , If gep ~ -E, on 
the other hand, the simple wave packet treatment is inaccurate |2^ . This fact is manifest also in our 
calculation: if ct£;p is not small enough, the integration in Eq. (|67D becomes problematic since Pic{E) 
is not necessarily meaningfully defined for negative E values. Anyway, there might be at least in 
principle a situation where the energy distribution of the neutrino wave function deviates considerably 
from a plane wave. Our calculation suggests that then the usual level crossing probabilities would 
not be totally reliable. 

Finally, it is to be noted that the use of some specific Pic{E) restricts the location of xq and x 
in Eqs. (Q) and (|62D; Plz{E), for example, is valid only if xo{x) is situated well before (after) the 
resonance. If, on the contrary, xq and/or x are/is in the resonance region, Plz{E) cannot be used. 
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5 Summary 



We applied the wave packet formalism in order to study neutrino oscillations in matter in the adia- 
batic limit, and found out that the effective oscillation and coherence lengths take into account the 
whole path the neutrino has traversed. Results for the linear and exponential density profiles were 
briefly presented. The corrections for the predictions of observable fluxes of solar neutrinos seem to 
be quite small. On the other hand, the matter may affect significantly the oscillation of the neutrinos 
e.g. in long baseline experiments for suitable values of parameters. 

We then considered the equation of motion for two neutrino flavors, and managed to solve it 
formally for arbitrary density profile. Our method clearly applies to any differential equation of the 
same form. 

Finally, we showed that the level crossing probabilities of the wave packets differ from those of 
plane waves. The difference is practically equivalent to a simple average over energy (cf. Eq. (|67D), 
not essentially related to the separation of wave packets. The finite width of the wave packet does 
not result in any observable effect for the physical situations we have considered. 

We could have continued our work by combining the wave packet treatment of Sec. 2 with the 
complete solution of the two flavor equation in Sec. 4.1. In that case the solution of Eq. (H) (for 
a = 1,2) should be replaced by Eq. (|^) in order to correctly take into account the effects due to 
nonadiabaticity. However, the presented results suggest that this would not reveal any new physics, 
so the calculation has been omitted. 

In this work we used a model that in principle is more physical and hence more accurate than the 
models used normally. On the other hand, since many of the calculations in this framework are quite 
complicated, we considered the limits of validity of the simpler plane wave approaches. Our results 
show that the present neutrino observations and the phenomena behind them can be described by a 
plane wave model accurately enough, regarding the current precision of the experiments. 
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Appendix A: Integration over x in Eq. ([T^) 

We can write Eq. (|T^ as 




(Al) 



where 




( 



/■^ dx' \ 

Jl Va{x') J 

j-F{x)), leading to 



(A2) 



The definition of the saddle point is in turn F'{xq) = {F'{x) 




Va{x') 



■^0 dx' 




xP 



L Va{x') 



■^0 dx' 



vl{L) = . 



(A3) 
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From this one can solve that 



va{x') alj^vm + alpvl{L) 

and 

-0 dx' _crM(0)(T-/o^^ 



(A4) 



(A5) 



and finally after some algebra 



1 (-^ -^0 vt{x 



2 



nxo) = -t / ,r 'r\. ■ (A6) 

Let us now examine the higher derivatives of F{x). The second derivative is 




In the relativistic limit Va{x) = 1 + 0{fil/E^), and Eq. (^) reduces to0 



v'Jx) . (A7) 



F"M = - ^ + O f ^;^^) . (A8) 

where effectively cr^p j^ = max{o"^p, o"^^}. From Eq. ( |A7| ) one sees also that the third derivative 
of F{x) includes terms proportional to v'^{x), {v'^{x)Y and v"{x). It is thus obvious that F"\x) 
and all the higher derivatives are negligible in the relativistic limit, and we can approximate (with 
F'(xo) = 0) 

F{x) ^ F{xo) + ^F"{xo) (x-xof . (A9) 



The relevant part of the integral of Eq. (|TJ) becomes 



J dxe^^""^ ^ e^(^o) J dxexp 



(AlO) 



where the final Gaussian integral gives only an unimportant numerical factor. 



^Thc adiabaticity further reinforces the smallness of v'J^x) 
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Appendix B: Additional calculations in connection with the 
two flavor equation of motion 

We present another way of solving Eq. Instead of Eq. (|5DD, we now define (compare to Ref. 



43 



^i(x) 



exp {t dx'^ 




exp (^-ij^^^dx'^ 



01 (x) 

02 (a;) 



leading to 



where 



'>2ix) 





iD*(x) 



-iD{x) 




(pi[x) 

02 (x) 



/ rx ,A(x')\ 

D{x) = 9'^{x) exp \-i J^^ ^^'^^j ■ 

Eq. (|B2]) yields easily 

0i(a;) = Ci(^l-jD-jDjD* + jDjD*jD + jDjD*jDjD* 

D I D* - 



+C2 1 + D 



D D* D 



D D* D D* + ■ 



02(x) 



Ci(^i + Jd*-Jd*Jd-Jd*JdJd* + J d* J d J d* J d 



-Co 1 



D* 



D* D + D* D D* + D* D D* D 



where the notation is as before (see Eq. (p3D). Finally one has 



exp t 



dx ^ 
xo AE 



V'i(xo) (i-JdJd* + JdJ D*J d j d* + - 



1p2ix) 



exp - 



D- 

A{x') 



D D* D 



'^li.Xo) 
+V'2(Xo) ( 1 



dx 

xo AE 
' D* 



D* D D* + 



D* / D 



D* D D* D 



where the energy dependence is not written down explicitly. We thus see that Eq. 
only form in which the solution of Eq. (^Up can be expressed. 

The equation of motion in the flavor basis is p4|, ^ (and many others) 



ipeix) \ _ _}_ ( -Aocos29 + Ac{x) Ao sin 261 
^^,{x) )~AE\ Ao sin 261 Aq cos26' - Ac(x) 



tpe{x) 
-^f^ix) 



(Bl) 

(B2) 
(B3) 



(B4) 



(B5) 
is not the 

(B6) 
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where Ac{x) = 2\/2EGFNe{x) and other notations are obvious. Since the calculation proceeds as 
above, we omit the details and content ourselves with giving the final answer: 



where 



i/jeix) = exp 



X 



ijj^ix) = exp 



X 



r dx'{AoCos20-A,{x')) 
M^o) (l- I F I F* + j F I F* I F J F* + 

^ r dx'{AoCos2e - A^{x')) 

^eK) (-^ J^*+'J ^* j ^ J ^* + - 



+V'^(xo) [l- J F* J F + J F* J F J F* J F + 

dx'iAo cos 26 - Ac{x')) 



Aosin2g 
F{x) ^ exp 



—I 
2E 



(B7) 
(B8) 
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